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Abstract 



, We show that factorization holds for the deeply virtual Compton scat- 

^SJ ! tering amplitude in QCD, up to power suppressed terms, to all orders in 

■ perturbation theory. The theorem applies to the production of off-shell pho- 
tons as well as real photons. We give a detailed treatment of the situation 

. where one of the two partons joining the parton density to the hard scattering 

Q,^ , has zero longitudinal momentum. 
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Typeset using REVT^ 



I. INTRODUCTION 



In this paper, we prove factorization for the deeply virtual Compton scattering (DVCS) 
amplitude in QCD up to power suppressed terms, to all orders in perturbation theory. This 
proof is important because of the recent great interest in DVCS |l|-[TO|. One important 
use of DVCS is as a probe of off-forward (or nondiagonal) distributions [lll,|3|,p!lHT3| . These 



differ from the usual parton distributions probed in inclusive reactions by having a non-zero 
momentum transfer between the proton in the initial and final state. 

A related process which is also used to probe off-diagonal parton densities is exclusive 
meson production in deep- inelastic scattering , for which a proof of factorization was 

given in |T2| . Compared with this process, DVCS is simpler because the composite meson in 
the final state is replaced by an elementary particle, the photon, and thus there is no meson 
wave function in the factorization formula. 

In the case of a scalar field theory, Anikin and Zavialov [|1^] proved a non-local operator 
product expansion, from which follows the factorization theorem for DVCS, as shown by 



Miiller et al. ^Tj. From the point of view of these papers, the new results in the present 
paper consist of an extension of the results to QCD. However, we do not derive an explicit 
form of an /^-operation for the coefficient functions, unlike Anikin and Zavialov. 

In a separate line of development, Ji [Q and Radyushkin provided key insights that 
indicate that a factorization theorem is valid for DVCS, and then Radyushkin provided 
an all-orders proof in 0. In this paper we provide an alternative proof, and give a new 
treatment of some problems that were treated in Ref. but that were perhaps not fully 
solved. The primary technical difference between Radyushkin's derivations and ours is that 
he uses the a-parametric representation for Feynman graphs, whereas we use the momentum 
representation, which we consider to be more direct. Our proof follows the general lines 
of proofs of factorization for other processes given in ||T8|JT9(] , and the most noteworthy 
feature is that, particularly for the case of production of off-shell photons, the proof is 
simpler than for any other process. Even for ordinary deep-inelastic scattering one needs 
to discuss the cancellation of soft gluon exchanges and of final-state interactions, whereas 
these complications are not present in the leading power for DVCS. 

The paper is organized in the following way: After stating the theorem in Sec. |T|, we show 



in Sec. Ill how the proof given in Ref. [12] is readily adapted to DVCS. The complications 



mentioned above concern the situation when one of the two lines connecting the parton 
density to the hard scattering carries zero longitudinal momentum, and these are given a 



detailed treatment in Sec. Ill H 



II. FACTORIZATION THEOREM 

The process under consideration is DVCS which is the elastic scattering of virtual pho- 
tons: 

^*{q) + P{p)^^^*\q') + P'{p-A) (1) 

where the diffracted proton P' may also be replaced by a low-mass excited state and the final- 
state photon can be either real or time-like. This process is the hadronic part of ep — > e'jp' 
for a real photon or of ep — > efi~^fi~p' for a time-like photon. 
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It is convenient to use light-cone coordinates with respect to the colhsion axi^. The 
momenta in the process then take the form: 




(2) 



Here, x is the Bjorken scahng variable, is the virtuality of the initial photon, is the 
proton mass, t = is the momentum transfer squared, and a is a parameter that specifies 
the virtuality of the outgoing photon: g'^ = aQ^. Thus, a = for a real photon and a > 
for a time-like photon. Finally, ~ means "equality up to power suppressed terms" . 
The theorem to be proved is that the amplitude for the process (j^) is: 



where fi/p is a nondiagonal parton distribution and Hi is the hard-scattering coefficient 
for scattering off a parton of type i. We let xi be the momentum fraction of parton i 
coming from the proton, so that X2 = Xi — x{l + a) is the momentum fraction which is 
returned to the proton by the other parton line joining the parton distribution and the 
hard scattering. There is implicit polarization dependence in the amplitude, n is the usual 
renormalization/factorization scale, which should be of order Q to allow calculations of the 
hard scattering coefficients within finite-order perturbation theory. The fi dependence of 
fi/p is given by equations of the DGLAP and Brodsky-Lepage kind [p^ JTT|JT^JT^JT3[] . The 
parton distributions in Eq. (|]), together with their evolution equations, are defined using 
the conventions of They may easily be transformed into those given in by a 

change of normalization and of kinematic variables. 




power-suppressed corrections, (3) 



III. PROOF OF THEOREM 



The proof of our theorem Eq. can be summarized as follows 0: 



We define a vector in light cone coordinates by: 




^ For a very detailed account of the basic steps and potential problems see Ref. |12|. 



3 




a) b) 

FIG. 1. a) Reduced graph for DVCS with direct couphng for the out-going photon to hard 
subgraph, b) The same without a direct hke couphng for the out-going photon. 

• Establish the non-ultra-violet regions in the space of loop momenta contributing to 
the amplitude. 

• Establish and prove a power counting formula for these regions. 

• Determine the leading regions of the amplitude. 

• Define the necessary subtractions in the amplitude to avoid double counting. 

• Taylor expand the amplitude to obtain a factorized form. 

• Show that the part containing the long-distance information can be expressed through 
matrix elements of renormalized, bi-local, gauge invariant operators of twist-2. 

A. Regions 

First let us establish the regions in the space of loop momenta contributing to the asymp- 
totics of the amplitude, i.e., the generalized reduced graphs. The steps leading to the gen- 
eralized reduced graphs are identical to the steps 1-3 in Sec. IV of Ref. [|l2l, i.e., scale all 
momenta by a factor Q/m, use the Coleman- Norton theorem to locate all pinch- singular 
surfaces in the space of loop momenta (in the zero- mass limit), and finally identify the 
relevant regions of integration as neighborhoods of these pinch-singular surfaces. 

The results are the two kinds of reduced graph shown in Fig. |l[ There, A and B denote 
collinear graphs with one large momentum component in the + and — direction respectively, 
H denotes the hard scattering graph, and 5* denotes a graph with all of its lines soft, i.e., 
in the center-of-mass frame all the components of the momenta in S are much smaller than 
Q. Note that, of the external momenta, p and p' belong to A, q' belongs to B oi H , and q 
belongs to H. 
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a) b) 

FIG. 2. Those reduced graphs that contribute to the leading regions in DVCS. 

When the two external photons have comparably large virtualities, the only reduced 
graphs are of the first kind, Fig. |I]a, where the out-going photon couples directly to the hard 
scattering. But when the out-going photon has much lower virtuality than the incoming 
photon, for example, when it is real, we can also have the second kind of reduced graph. 
Fig. |l|b, where the out-going photon couples to a i? subgraph. As we will see later, power 
counting will show that the second kind of reduced graph. Fig. ^d, is power suppressed 
compared to the first kind, with a direct photon coupling. This implies that we will avoid 
all the complications which were encountered in |12[ that are associated with the meson 
wave function. 



B. Power Counting 

Each reduced graph codes a region of loop-momentum space, a neighborhood of the 
surface vr of a pinch singularity in the zero-mass limit. The contribution to the amplitude 
from a neighborhood of vr behaves like Qp^'^\ modulo logarithms, in the large-Q limit, with 
the power given by 

p(7r) = 4 — n{H) — ^^(quarks from S to A, B) — 3#(quarks from S to H) 

-2#(gluons from S to H). (4) 

where n{H) is the number of collinear quarks, transversely polarized gluons, and external 
photons attaching to the hard subgraph H. Such results were obtained by Libby and Sterman 
P(]| , PT| . The particular form of Eq. was given in |^ together with a proof that applies 
without change to DVCS. 



The well-known problem of gluons with scalar polarization (see, for example, |T8tp2| , p3[ ) 
will be dealt with later on. Suffice it to say here that gluons with such a polarization can 
be factorized into the parton distributions by using gauge-invariance arguments. 
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FIG. 3. Soft gluon loop attaching to collinear line. 
C. Leading Regions 

The leading regions for the amplitude are those with the largest exponent p(vr) in Eq. 
(|). It is easy to see that these correspond to the reduced graphs in Fig. ^ independently of 
whether the out-going photon is real or far off-shell. The corresponding power is Q'^. These 
reduced graphs have direct photon couplings to the hard subgraph, they have exactly two 
parton lines connecting the collinear subgraph A to the hard subgraph H, and they have no 
soft lines connecting to H. The two kinds of graph differ only by the absence or presence of 
a soft subgraph that connects to A alone. 

Among the other reduced graphs, which are non-leading for our process, are those of the 
type in Fig. ^b, which are leading in the case of diffractive meson production, where the 
leading region gives Q~^. 

In the case of a photon that is off-shell by order Q^, the amplitude for production of the 
photon behaves like Q^, the same as for a real photon. However, the physically observed 
process includes the decay of the time-like photon (to a for example), which results in 

a power suppression of the observed cross section by l/Q"^ compared with the cross section 
for making real photons. 



D. Proof of absence of a soft part in leading regions 

As mentioned in Sec. [Ill C[ there might, in principle, be a soft part S in the leading 
reduced graph connected solely to the A graph by gluons, as shown in Fig. Note that 
by Eq. (^), quarks connecting S to A would lead to a power suppression. We will now show 
that this soft part S is indeed absent, and so we only need to consider regions of the form 
of Fig. |a. 

We will first examine a simple one loop example. Fig. ^. The external quark is part of 
the A subgraph in Fig. |^, and the gluon is soft. So we parameterize the momenta by: 

kA = {xiP^,kJ,kA,i.), 
k = {k+,k-,k±). (5) 

where the fcj is 0{Q) and all the other components are of 0{m) or smaller. 

If we omit irrelevant factors in the numerator, the loop integral takes the following form: 



soft k {k? + ie) [{kA — ky — m? + it] 
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d^k- 



1 



soft k {2k+k- -k\ + ie) 2{xip+ - k+){kj - A;") - - ki_f - m? + ie 

1 



dkA 



soft k {2k+k- -k\ + ie) 2xip+{kj - k^) - {kA,L - k±y - w? + ie 



(6) 



As before, ~ stands for "equality up to power suppressed terms". As one can see, there is 
no fc^-pole in the second part of the denominator and we can freely deform the contour in 
k'^ to avoid the pole in the soft gluon propagator. This takes us out of the soft region for k. 

In the general situation. Fig. we can use a version of the arguments in Ref. |I2|J2^ 
to show that the soft momenta kf can be rerouted in such a way as to exhibit a lack of a 
pinch singularity. The essential idea is that one can find a path backwards or forward from 
one external line of S to another external line of 5*. The loop is completed along lines of A, 
all of which have much larger + momenta than what is typical of soft momenta, and hence 
there is no pinch. 



E. Subtractions 



The subtractions necessary to avoid double counting in the amplitude are treated exactly 
the same fashion as the ones in Sec. VI of Ref. |jl2|, since the distributional arguments to 



construct the subtraction terms on a pinch-singular surface vr presented there are very general 
in nature and are not limited to the case of diffractive vector meson production that was 
considered in 



12 



The above statement leads to the following asymptotic form of the amplitude 

AsyT = SrAsyT = Ax H. 
where T stands for a possible graph for the amplitude T. 



(7) 



F. Taylor expansion 

We now obtain the leading term in the hard subgraph, when it is expanded in powers of 
the small momenta. The arguments used are exactly analogous to the ones used in Sec. VII 
of Ref. [|12| except that we do not have to deal with a B subgraph as was the case in [jl2 



So we have: 

A X if ~ j dk\H(q, q\ {k+, 0, 0^), (A+ - k+, 0, 0_ 

I dkjd^kA,±A{kA,A-kA), (8) 

where fc^ is the loop momentum joining the A and H subgraphs, and again ~ means equality 
up to power-suppressed corrections. Eq. has already a factorized form. However we still 
have to deal with the extra scalar gluons that may be exchanged between the subgraphs A 
and H; this will be done in the next subsection. 
Eq. M) can be written in the following way: 
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a) 



b) 



FIG. 4. a) A scalar gluon attaching the colhnear subgraph to the hard subgraph H in the 
unfactorized form, b) Factorized form after application of gauge invariance and Ward-identities. 
The double line represents the eikonal line to which the scalar gluon attaches. 

AxHc^T.,j dk+Qiq,q',k+)0,{p,p',k+), (9) 

where the Ci are the short distance coefficient functions and the Oi are the matrix elements 
of renormalized light-cone operators. 



G. Gauge Invariance 



In order to identify the Oj with the parton distributions as defined in |]T2| (for example), 
it is necessary to show that all gluons with scalar polarization attaching to the hard graph 
can be combined into a path-ordered exponential. Fig. |^ shows the example of one scalar 
gluon. This was shown in 

In this way 



Ref. m 



|P3[]. Another version of a proof was given in Sec. VII. D of 
which relies on very general results recently obtained by Collins 



we obtain exactly the same parton distributions as in [|T2|, namely: 



fq/p 



oo 



dy- 



-IX2P 



OO An 

°° dy 



'^"(p|TV^(0,y-,Ox)7+^V^(0)b'), 



— oo 



27r XiX2P'^ 



'tX2P^y 



(p|TG+(0,2/-,Ox)W^+(0)b'). 



(10) 



Here, V represents a path-ordered exponential of the gluon field that makes the operators 
gauge invariant. The variable X2 is the same as in Eq. (^. The evolution equations are the 
same as in [Q,R1 and [IT^l. 



H. Partons with = 0: breakpoints and endpoints 

In the factorization theorem Eq. (0), the integral over the fractional momenta includes 
the points Xi = and X2 = 0. At these points, the hard scattering coefficient for DVCS 
has a pole, and so we appear to get a logarithmic contribution to the cross section from a 
region in which one of the lines joining the parton density to the hard scattering subgraphs 
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FIG. 5. Particular example of potentially problematic diagram. 



is soft instead of coUinear. This apparently contradicts our power-counting result that such 
a region gives a non-leading power. This phenomenon was investigated by Radyushkin 0. 
In this section, we will use momentum-space methods to give a general demonstration that 
the region in question does not give a problem. 

First, let us observe that the region of integration over xi in the factorization formula 
Eq. (^) is — 1 + x < Xi < 1. This is proved by the methods of light-front perturbation theory, 
by requiring that the intermediate states in Fig. ^ be physically allowed. See Ref. PJTB[] for 
detailed derivations and discussions. The points Xi = and X2 = at which the potential 
problem arises are what we will call "breakpoints", since they occur in the middle of the 
range of integration where one of the two lines changes direction. 

We continue by examining a particular case, illustrated in Fig. |^, and showing how the 
argument generalizes. To simplify the example, let us restrict our attention to regions where 
the subgraph A and the lower three lines (p, p' and k — p') have their momenta collinear to 
the proton. We will also require the two quark lines, k + p — p' and k, on the sides of the 
ladder to have their momenta either collinear to the proton or soft. The example is very 
similar to one treated by Radyushkin in Ref. |0. 

We will also only need the case of the production of a real photon, g'^ = 0, since this is 
where the problem arises. 

The top loop of the graph has the form: 

/■ ,4, Numerator factors 

J [k? — m? + ie) [{k — p'Y — m? + ie] [{k + p — p'Y — rn? + ie\ [{k + q'Y — rn? + ie\ 

(11) 

When both k and k + p — p' are collinear to A, the top line is off-shell by 0{Q'^), and it is 
correct to use the collinear approximation 

(12) 



{k + q')"^ — m"^ + ie X2Q'^/x + ie' 



where X2 = k'^ /p^- A corresponding replacement is also to be made in the numerator 
in U . The right-hand side of Eq. (12) exhibits the afore-mentioned pole at X2 = 0. The 



result of applying the coUinear approximation is to give the appropriate contribution to the 
factorization formula Eq. (|]). 

The coUinear approximation becomes invalid when k becomes soft, i.e., when X2 0. 
We must now demonstrate two facts. The first is that, when is in a neighborhood of the 
soft region, the coUinear approximation is valid after integration over k. The second fact is 
that the use of the coUinear approximation does not give an important contribution from 
some other region of k which is absent in the original, unapproximated graph. We now 
examine the integral U in the neighborhood of the soft region for k. It has the following 
form: 

r 4 Numerator factors , , 

' " Xoft k [2k+k- -kl-m^ + te] [2p'+{p'- - k-) - {p' - k)l - w? + le] ^ ^ 

1 

[2{p+ - p''^){p~ - p'^ + k-) - {p-p' + k)l_ -m? + ie] -k\-m? + ie) ' 

where we have neglected k~^ in the collinear-to-A lines, k~ in the collinear-to-i? line and 
X ^ Xi- We have also ignored the numerator factors, which are an irrelevant complication 
for our purposes. 



According to the power-counting results of Sec. |1IIB| , which are obtained from [£4[, the 
soft region for k gives a power-suppressed contribution. This estimate assumes that all 
components of k are comparable (in the Breit frame), and is obtained as follows. Let the 
magnitude of the components of k be m. Then the order of magnitude of the soft part 
of f/ is a product of factors l/{nv'Q^) from the denominators, from the phase space, 
and Q'^m? from the numerator, for an overall power m/Q. This result can be obtained by 
writing down the largest components in the trace and propagators of Fig. ^ and Eq. (|13|). 
Moreover in this region it is correct to replace the fourth propagator in Eq. (|l^) by its 
coUinear approximation Eq. (|1^), so that we do not lose the factorization theorem. 

However, if the components of k are asymmetric, this estimate no longer holds. In par- 
ticular if the longitudinal components of k are of order k^ ~ /Q while the transverse 
components remain of order m, then we get contributions of order 1/m^ from the denomi- 
nators, m^/Q^ from the phase space, and Q'^m? from the numerators, for a total of m^Q^. 
This shows that the contribution from this region is unsuppressed for large Q. 

At this point we must appeal to the contour deformation arguments of Ref. P3. It is 



only when the integration over k is pinched in the region in question that it needs to be taken 
into account. In the dangerous region we have k~^k~ k\, so that the only k~^ dependence 



in Eq. (13) is the pole in the fourth denominator. We can therefore deform k~^ into the 
complex plain a long way out of the region we are considering, indeed all the way to the 
collinear-to-A region. Then the collinear approximation is valid so that we can replace the 
graph by its contribution to the factorization formula. 

This contour deformation argument is completely general, as explained in Sec. HIE of 
Ref. [24|. Whenever we have a soft momentum with k~^k~ -C k\, the contour of k~^ can be 



deformed away from poles in the jet subgraph associated with with the produced photon. 
Since all the relevant singularities are in the final state, they are all on the same side of the 
real axis. 
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Now that we have estabhshed in more detail that the only leading regions are those 
symbolized in Fig. §(a), we can apply the collinear approximation as described earlier, and 
hence we obtain the factorization theorem. 

But we still see the following problem. In the factorization theorem, Eq. @, the parton 
densities are non-analytic at the breakpoints xi = and X2 = 0, whereas the coefficient 
function has a pole at each of these points. Again consider the collinear approximation to 
Fig. 1^ in the region we were considering. The parton density is non-analytic when X2 = 0, 
while the coefficient function has a pole there, as is seen from the right- hand- side of Eq. 
(|I^). So we cannot literally apply the contour deformation argument. 

What we will show is that the parton density is continuou^ at the breakpoint, so that 
it can be written as the sum of a function that is analytic at X2 = and a function that 
has a zero at X2 = 0. The only potential leading twist contribution near the breakpoint 
is associated with the non-zero analytic term to which the contour deformation argument 
applies. 

To prove this property of the parton density at a breakpoint, consider a general graph 
for the parton density, as shown in Fig. |^. We have found it convenient to change the 
labeling of the momentum compared with the previous figure. As always, the k~ and 
components of k have been short circuited and are integrated over. The fc-line gives a pole 
at k~ = {k'j_ + — ie)/2k~^ = {k\ + m? — ie)j2x\P^^ while the k + q' — g-line gives a pole at 
k- = [{k + q- + m?- ie]/2(A;+ - ^P+) = [{k + q- q')\ + - ie]/2a;2P+. Here, ^ is the 
fractional longitudinal momentum transfer 1 — P'^ / . In addition there are poles from the 
collinear-to-y4 in the blob. For example if the blob consists of a single line, we have a pole at 
k- = p--{P+ml-ie)/2{l-xi)P+ or at fc" = -P'~ + [{k+P'f+ml-ie)/2{l-^+xi)P+. 

As we vary xi, the k~ contour can generally be deformed to avoid the poles, so that 
we have analytic dependence on Xi. The possible exceptions occur when the k~ contour is 
pinched for finite k~ or when a singularity coincides with the endpoint of the integration 
at k~ = oo. A pinch never occurs; in the general case this is a consequence of the Landau 
rules. But endpoint singularities occur, and these are precisely at the breakpoints. 

For example if fc"*" — > 0, then fc"*" can approach from above and below. The pole giving 
us trouble stems from the fc-line, all other propagators are unproblematic in this case, since 
their poles are at finite k~ . The pole in k~ approaches +oo — ie as one approaches 0+ and 
—oo + ie as one approaches 0~. This means that the k~ pole crosses the real axis at infinity. 
Hence the parton distribution is non-analytic there. Since the singularity is at \k~\ = oo the 
other propagators have large denominators, and hence we get a zero for the non-analytic 
part of integral at the breakpoint. Thus the parton density is continuous at the breakpoint, 
as claimed. This result enables the factorization formula to be valid in the neighborhoods 
of the breakpoints. Since the other poles in the k~ integral are on opposite sides of the real 
axis, the parton distribution is non-zero at the breakpoints. 

Effectively the crossover of the pole occurs when is in a coUinear-to-i? region, which we 
know is power suppressed. This indicates that the argument we have just given generalizes 
to all graphs. 



^ Radyushkin indicated briefly in Ref. how such a property is to be proved from the a repre- 
sentation for his double distributions. 
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k+q-q' 




FIG. 6. Parton distribution amplitude. 



We also remark on the behavior at the endpoints. Let us look at the case fc"*" p'^. We 
find that another of the poles pole runs off to —oo this time and crosses the real axis there. 
But now all the other poles are on a single side of the real axis, so that the sole contribution 
to the parton density comes from the pole at infinity, and hence there is a zero of the parton 
density at the endpoint. 



I. Completion of Proof 

Using the definitions of the parton distributions and the hard scattering coefficients we 
finally obtain Eq. Note that the theorem is valid for the production of a real photon 
which directly goes into the final state and for the production of a time-like photons that 
decays into a lepton pair. 



IV. CONCLUSION 



We have proved the factorization theorem for deeply virtual Compton scattering up to 
power suppressed terms to all orders in perturbation theory. The form of the theorem is 
independent of the virtuality of the produced photon. 

After this work was complete, we learned from X.-D. Ji that he and Osborne have also 
constructed a proof of factorization for DVCS PB . 
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